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摘 要 :采用 精确 解析 求解 方法 研究 石墨 烯 片 (GPL) 增 强 功 能 梯度 梁 的 非 线 性 屈曲 及 后 屈曲 力学 
行为 。 研 究 中 假定 GPL 均匀 地 分 布 在 梁 的 每 一 层 , 但 其 重量 分 数 沿 厚度 方向 呈 梯 度 或 均匀 变化 , 导 
致 结构 的 物性 参数 沿 厚 度 方 向 连续 梯度 变化 ,并 根据 Halpin-Tsai 力学 模型 和 混合 率 法 则 得 到 材料 
,的 等 效 物性 参数 。 同 时 基于 Euler-Bernoulli 梁 理 论 和 Von-Karman 非 线 性 应 变 位 移 关 系 由 能 量变 分 
法 获得 石墨 烯 片 增强 功能 梯度 梁 的 非 线 性 控制 方程 ,对 此 方程 进行 精确 解析 求解 获得 临界 载荷 、 慑 
曲 模 态 以 及 后 属 曲 变形 的 解析 表达 式 。 并 进行 参数 研究 分 析 了 GPL 的 重量 分 数 、 几 何 形状 及 尺寸 
和 分 布 方 式 对 该 复合 材料 梁 非 线 性 屈曲 和 后 屈曲 的 影响 。 
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Exact solutions for buckling and post-buckling of functionally 
grade graphene reinforced composite beams 
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2. Liuhuan Transmission ( Xian) Technology Co.,Ltd,710000 Xi’an, China) 


Abstract : In this paper, nonlinear buckling and post-buckling behaviors of functionally graded beams rein- 
forced by graphene platelets (GPL) are studied by exact analytical method. It is assumed that GPL is uni- 
formly distributed in each layer of the beam, and its weight fraction varies gradient along the thickness di- 
rection , which leads to continuous and gradient changes of the materials properties along the thickness. Ac- 
cording to Halpin-Tsai mechanical model and mixing rate rule ,the equivalent materials properties are ob- 
tained. Simultaneously , based on Euler Bernoulli beam theory and Von-Kérman nonlinear strain-displace- 
ment relationship, nonlinear governing equations of the graphene platelets reinforced functionally graded 
beam are obtained by energy variational method. And the analytical expressions of critical buckling loads, 
buckling modes and post-buckling deformations are obtained through exactly and analytically solve the 


nonlinear governing differential equations. The effects of the weight fraction, geometrical parameters and 
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distributions of GPL on the nonlinear buckling and post-buckling of the composite beam are analyzed by 


parameters research. 
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Ay ois his TURE RI Ea PE oe BY JI A PE 
能 ,使 它 成 为 非常 理想 的 增强 材料 。 将 少量 的 石墨 
烯 或 其 衍生 物 作为 增强 体 加 入 基体 成 为 石墨 烯 增强 
复合 材料 , 比 传统 碳纤维 和 碳 纳米 管 更 能 显著 提 
高 结构 的 强度 和 刚度 ”1 。 将 功能 梯度 材料 ”1 卓越 
的 设计 理念 推广 到 石墨 烯 增强 复合 材料 中 ,也 就 是 
K-A Sei Fr (graphene platelets, GPL) 县 成 的 小 块 作 
为 增强 体 方向 随机 且 均 匀 地 散布 在 每 层 基 体 中 , 石 
墨 烯 含量 在 基体 厚度 方向 上 逐 层 呈 梯度 或 均匀 分 
布 ,形成 新 一 代 的 先进 复合 材料 , 即 功 能 梯度 石墨 烯 
增强 复合 材 Æ} U (functionally graded graphene plate- 
loreinforeed composite, FG-GPLRC ), FG-GPLRC 
也 县 有 非凡 的 力学 性 能 , 主要 取决 于 石墨 烯 , 可 以 通 
壕 会 理 设计 石墨 烯 片 的 分 布 方式 .重量 分 数 和 几何 
参数 得 到 最 优 的 整体 力学 性 能 。 

QOFC-GPLRC 结构 的 力学 性 能 研究 近年 来 已 成 为 
国体 力学 中 一 个 新 的 研究 领域 。 在 线性 稳定 性 方 
HL KITIPORNCHAI a RH EZX GPL 增强 功 
能 梯度 多 孔 梁 的 自由 振动 和 弹性 届 曲 进行 了 研究 。 
WO 等 5 对 处 于 热 环境 中 同时 受 轴 向 载荷 作用 的 多 
层 在 墨 烯 增强 功能 梯度 梁 的 动力 稳定 性 进行 了 研 
FE YANG 等 9 基于 能 量 法 研究 了 石墨 烯 增强 纳米 
复 吹 材料 拱 在 点 载荷 作用 下 的 动态 届 曲 。LU 等 ” 
基于 一 阶 剪 切 变形 板 理论 和 最 小 势能 原理 研究 了 石 
墨 医 增强 功能 梯度 压 电 板 的 届 曲 特性 。KIANIW 对 
热 环 境 中 受 外 部 压力 作用 下 功能 梯度 石墨 烯 增强 锥 
壳 的 届 曲 进行 了 研究 。TAM 等 小 对 含有 开放 边缘 
裂纹 的 石墨 烯 增强 功能 梯度 复合 材料 梁 的 屈曲 进行 
了 研究 。SHAHGHOLIAN 等 5 基于 一 阶 剪 切 变形 
理论 对 石墨 烯 增强 多 孔 纳米 复合 材料 圆柱 过 的 屈曲 
进行 了 研究 。PARASHAR 等 '" 采用 有 限 元 法 对 受 
轴 向 载荷 作用 下 的 石墨 烯 增强 环 氧 复合 板 的 屈曲 进 
行 了 研究 。 通 过 研究 表明 当 石 墨 烯 片 的 体积 分 数 是 
6% 时 ,该 结构 的 屈曲 强度 提高 26%。EBRAHIMI 
等 5 研究 了 热 环 境 中 GPL 增强 聚合 物 复 合 材料 板 
在 黏 弹性 地 基 上 的 屈曲 特性 。 同 时 考虑 了 GPL 的 4 
种 分 布 模式 ,通过 Hamilton 原理 建立 了 GPL 增强 板 
的 控制 方程 ,用 Navier 解法 对 该 问题 进行 了 分 析 求 
fiko SUN 等 基于 Donnell 壳 理 论 对 石墨 烯 增强 功 
能 梯度 复合 材料 圆柱 壳 的 扭转 屈曲 进行 了 研究 ,并 


由 此 得 到 了 FG-GPLRC 壳 体 的 临界 屈曲 载荷 和 相应 
的 屈曲 模 态 。 

以 上 研究 仅 限于 线性 问题 。 由 于 FG-GPLRC 结 
构 的 高 强度 .高 韧性 使 其 能 承受 较 大 变形 ,对 该 类 结 
构 进行 非 线性 大 变形 分 析 远 比 线性 小 变形 分 析 更 加 
有 意义 。 在 非 线性 稳定 性 的 研究 方面 ,YANG 等 中 
在 考虑 弹性 地 基 影 响 的 同时 通过 微分 求 积 法 对 石墨 
烯 增强 功能 梯度 复合 材料 梁 的 屈曲 及 后 屈曲 问题 进 
行 了 研究 。CHEN 4U php Fy ELBE RE ARIE OTE T 
沫 金属 基 GPL 增强 多 层 功 能 梯度 多 孔 纳米 复合 材 
料 梁 的 后 屈曲 行为 。SHEN 等 "研究 了 功能 梯度 石 
墨 烯 增强 复合 材料 层 合 梁 在 弹性 地 基 上 的 热 后 屈曲 
及 非 线 性 弯曲 。WANG 26") 对 GPL 增强 压 电 复合 
材料 梁 的 屈曲 和 后 屈曲 行为 进行 了 分 析 , 考虑 剪 切 
变形 和 非 线性 应 变 -位 移 关 系 ,利用 虚 功 原理 建立 了 
该 梁 的 控制 微分 方程 ,采用 微分 求 积 法 定量 地 研究 
影响 届 曲 性 能 的 因素 。ANSARI 等 :的 基于 一 阶 剪 
切 变形 过 理论 和 Von-Karman 非 线性 几何 关系 研究 
了 功能 梯度 石墨 烯 片 增强 聚合 物 纳米 复合 材料 圆柱 
壳 的 后 届 曲 ,并 利用 Fourier 级 数 和 微分 求 积 数 法 提 
出 了 半 解 析 解 的 方法 。SONG 等 i' 在 考虑 弹性 地 基 
和 边缘 裂纹 影响 的 同时 基于 一 阶 剪 切 变形 梁 理 论 ， 
研究 了 石墨 烯 复合 材料 多 层 梁 的 热 届 曲 和 后 屈曲 。 
SAHMANI £% 基于 非 局 部 应 变 梯度 理论 研究 了 轴 
向 载荷 作用 下 GPL 增强 多 层 功 能 梯度 复合 纳米 板 
的 非 线性 屈曲 及 后 屈曲 。 

目前 石墨 烯 增强 功能 梯度 材料 结构 非 线 性 力学 
行为 的 已 有 研究 还 比较 少 ,尤其 对 该 结构 非 线性 稳 
定性 问题 的 非 线性 控制 微分 方程 进行 精确 解析 求解 
方面 鲜 有 文献 报道 。 本 研究 在 考虑 非 线性 大 变形 的 
情况 下 ,基于 Halpin-Tsai 力学 模型 2 ,用 解析 法 精 
确 求解 石墨 烯 增强 功能 梯度 粱 的 非 线性 稳定 性 控制 
方程 ,获得 临界 屈曲 载荷 . 届 曲 模 态 和 后 屈曲 变形 的 
表达 式 ,并 分 析 GPL 的 重量 分 数 .不 同 的 分 布 模式 、 
几何 形状 和 尺寸 等 对 临界 屈曲 载荷 的 影响 。 


1 数学 模型 


1.1 材料 物性 参数 
考虑 图 1 ERKA LIEN b AN h 的 石墨 烯 增 
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成 ,并 且 每 层 中 的 石墨 烯 为 均匀 分 布 ,但 不 同 层 中 石 
墨 炳 含量 不 同 。 本 研究 考虑 各 层 中 石墨 烦 的 含量 具 
有 U-GPLRC ,X-GPLRC 和 0-GPLRC 3 种 不 同 分 布 
模式 如 图 1 所 示 , 图 中 颜色 越 深 则 石墨 烯 浓度 越 大 。 
U-GPLRC 表示 GPL 在 每 一 层 均 匀 分 布 且 浓 度 相同 ; 
X-GPLRC 表示 GPL 在 上 层 和 下 层 具有 较 高 的 浓度 
而 在 中 间 浓 度 最 低 ;0-GPLRC 表示 GPL 在 上 、 下 层 
有 具 有 较 低 的 浓度 而 中 间 层 具有 较 高 的 浓度 。 这 3 种 
分 布 模式 的 第 下层 GPL 体积 分 数 Vem 表达 分 别 为 
U-GPLRC: Vin, = Veo. » 

X-GPLRC; Vi, = 2V [2k —N, -1|/NL. 
O-GPLRC; Von, = Vén - |2k - N, — 197M. 


IF.: k = 1,2,3…N ，N 是 梁 的 总 层 数 ; Voy, 是 
总 的 体积 分 数 , 可 表示 为 
© p Sori 
© Vopr = 
© gcpL + Cpcp./pm? 1 = gcpL) 
HD: 8 cpL 是 石墨 烯 的 质量 分 数 ; Pu All Pop. 分 别 为 
人 墨 烯 的 密度 。 


U-GPLRC X-GPLRC 


O-GPLRC 


图 1 3 种 不 同 的 GPL 分 布 模式 
Fig.1 Three different GPL distribution patterns 


FG-GPLRC 结构 第 层 的 弹性 模 量 £5” 基于 目 
前 广泛 采用 的 Halpin-Tsai 微观 力学 模型 给 出 "3 , 即 


EP = ( 31 + ENVE 51 Ua 
c cin + Ck M 
8 1- NLV cr 8 1- Mw Verr 


chinaXiv:: 


(1) 
> (Eom/Ew -1 _ CE op/Ey? ~ 1 
oe m= CEcp /Ey) + éL’ AS CE op /Ey) + Ew : 
Ey 和 E cr 分 别 为 基体 和 石墨 烯 的 弹性 模 量 ， ÉL All 
Ey 表示 GPL 的 几何 形状 尺寸 参数 ,具体 为 
= lop. _ af Scr. 
é = (2), by 2( a) 
hop 和 bcp 分 别 表 示 GPL 的 长 度 宽度 和 


hep. 

式 中 /cr 
厚度 。 

FG-GPLRC 的 泊 松 比 可 由 混合 定律 给 出 


J 
a fE m hua mR NAO ERT 591 


Ck) ck k 
Vo = Vor Ver, + UML — Vege) (2) 


式 中 Vy 和 VepL 分 别 为 基体 材料 和 石墨 烯 的 泊 松 比 。 
12 梁 的 控制 方程 


根据 Euler-Bernoulli 粱 理论 , 同时 考虑 Von- 
Kármán 非 线性 应 变 位 移 关系 ,可 以 得 到 梁 的 轴 向 应 
变 表达 式 为 
e, = +t (%) + 2K (3) 
式 中 ; x 和 = 分 别 为 轴 向 和 横向 坐标 ; w A w 分 别 为 
铀 向 位 移 和 找 度 ; x = - CF 为 梁 中 性 层 的 曲率 。 


考虑 线 弹 性 范围 内 的 变形 ,由 弹性 本 构 关 系 可 
1345 k JÆ FG-GPLRC WIEN A! 
Ch _ Ec » du 1 dw 2 dw 
| 
将 其 沿 厚 度 方 向 进行 积分 ,可 获得 粱 横 截 面 上 
的 轴 力 N, MEE M, 表达 式 , 即 


N, = foda = 4, [42+ (4) ] Ga) 
M sf dA =- D, dw 
= Fae =-= DV die (Sb) 


式 中 : 4 为 梁 横 截面 面积 ; 4， 和 Di 分 别 表示 梁 的 
拉 压 刚度 和 老 曲 刚度 ,通过 积分 计算 可 得 


N, ck 
Ec h po 
{A,, Dy} of mi :| 
iiaii pa Lao? N? 3 


采用 最 小 势能 原理 67 = 8U + V) = 0 可 以 推 
导出 梁 的 控制 方程 , 式 中 开 、V 和 了 分 别 表示 系统 
的 总 势能 ,应变 能 和 外 力 势 能 , 8 表示 变 分 。 

FG-GPLRC 梁 的 应 变 能 及 其 一 阶 变 分 为 


1 1y i 
U = ratat = 和 Ik N 8e? + M.8x)dx 


spet = Sg F(t)” 。 考 虑 梁 上 作用 给 向 压缩 


V = - [4p( 22) ax see P Pes( te) ax 
小 势能 原理 得 到 梁 的 控制 方程 为 


dN, 
ae = 0 (6) 


ae 
(et 
地 


dw Ëw _ 
-Dı (fee 0 (7) 
对 FG-GPLRC 深 的 屈曲 及 过 屈曲 问题 的 研究 ， 
即 对 方程 组 (6) 和 (7) 进 行 求解 ,还 需 同 时 考虑 边界 
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条 件 。 在 本 研究 中 考虑 梁 的 两 端 均 是 夹 紧 边 界 , 相 
应 的 边界 条 件 为 


00 (8) 


1.3 精确 解析 求解 


对 给 定 边 界 条 件 的 非 线性 控制 微分 方程 边 值 问 
题 进行 精确 解析 求解 。 首 先 将 式 (6)、 式 (7) 采 用 位 
移 形式 表示 为 


ll (ls w 
人 
(10) 
二 再 将 式 (9) 对 x 进行 一 次 积分 可 得 
fe, A, [4 + 3 (%) | = B, (11) 
加 二 次 积分 为 
pa Anfu + [> (e a] = Bx+B, (12) 


OTR B, M By ,将 边界 条 件 (8) 代 人 式 
(12) 可 得 
-一 再 将 B, 所 得 值 代入 式 (11 ) 得 到 
~ du 1 dw 2 E An 1 dw z 
Saltz la) I= (we 
-三 将 上 式 整 体 代入 式 (10) 并 整理 即 可 得 梁 挠 度 
的 控制 方程 , 即 
dfw > Anr dw\? dw = 
rd ae) tla a 
为 了 使 所 求 的 解 具有 一 般 性 ,引入 无 量 纲 量 , 即 
w Anh? PP Dy 
a aa, T ， “Da ° Da? 


通过 代入 无 量 纲 量 , 可 获得 石墨 烯 增 强 功 能 梯 
度 梁 拨 度 的 无 量 纲 控制 方程 , 即 
dw F rydW gd 
+ [P oc (Sy) dx] = 9 (14) 


dX dX’ 
夹 紧 端的 无 量 纲 形式 边界 条 件 为 
dW 
X=0,1h,W = 7 =0 


为 了 使 所 求 的 方程 形式 上 更 加 简洁 , 令 


parE fafa a 
则 式 (14) 转 化 为 
dd -0 (16) 
方程 (16) 的 通 解 为 
W(X) = Cicos(AX)+ Cysin(,X) + CX + C, 
(17) 


AP C, C, CG 和 C, 表示 积分 常数 。 进 一 步 将 上 
述 边界 条 件 代 入 式 (17) 中 有 
Wo) = C, + C, = 0 
W(1) = C,cosA + CysinA + C, +C, = 0 
(18b) 
(18c) 


(18a) 


WO)= CA+C =0 
W (1) =- CAsinA + C,AcosA + C; = 0 


(18d) 

令 C, = C , 联 立 求解 式 (18a) ~ (18d) 可 得 

B cosà - 1 _ (1 - cosà)À 
r= a Os ee sind 一 从 fks 


此 时 ,精确 解析 求解 获得 了 在 该 边界 条 件 下 石 
墨 烯 增强 功能 梯度 梁 届 曲 问 题解 的 形式 , 即 
cosA - 1 


sind 一 从 


WX) = cl — cos( AX) + sin (AX + 


ACE = eh) 
WP C 为 与 轴 辐 载荷 已 有 关 的 常数 。 再 将 式 (19 ) 
代入 式 (15 ) 可 得 
pan ee (20) 
JF NA 
替换 式 (18d) 中 的 Ci, Ca AC, 可 得 边界 条 件 
为 两 端 夹 紧 时 ,石墨 烯 增强 功能 梯度 梁 届 曲 的 特征 
方程 , 即 


2 -2coshA -AsinA = 0 (21) 

由 特征 方程 (21) 可 以 求 得 边界 条 件 为 两 端 夹 
紧 时 的 无 量 纲 特 征 值 为 

A, = 6.283 , A, = 8.987 ,A, = 12.566, A, = 
15.451 , A; = 18.850 ,--- 

采用 归 一 化 方法 令 C = 1 KRME A, RAR 
(19) 可 得 石墨 烯 增强 功能 梯度 梁 的 第 n 阶 届 曲 
模 态 。 
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(cosA, — 1) CsinA,X - AX) 
sind, — A 


W, = (1 - cosd,X) + 


(22) 
由 最 小 特征 值 A，= 6.283 , 即 为 27 ,将 其 代入 
到 式 (15 ) ,可 得 两 端 夹 紧 时 FG-GPLRC 梁 的 屈曲 临 
界 载荷 的 解析 表达 式 , 即 
D 


P, = 4a — 
cr T Do 


tif Dy = [一 Ya de Jue ASS A HUE 


对 齐 次 方程 (16) UE TOR AT Se 
Fc-CPLRC 染 的 临界 届 曲 载荷 对 应 的 后 届 曲 挠 度 表 
达 式 为 


X) =— = 5 天 -10 = cos(27X) 1 
(23) 
一 数值 算 例 分 析 
Siti clades anes, ee 


Seng ae ton 
OD TPa 密度 为 1.062 g/ em’ 、 泊 松 比 为 0.186。 如 
FREI BLO ,石墨 烯 的 质量 分 数 为 0.3% ,几何 尺寸 为 
ee 25 em bom =15hmn 和 jn = 1.5 pm ™!, 

”为 了 证 明 本 研究 的 理论 推导 及 求解 过 程 的 正确 
性 ,将 本 研究 的 结果 与 文献 [14] 中 使 用 微分 求 积 法 
所 得 的 相应 结果 作 比 较 。 引 入 文献 [14] 中 的 无 量 


Shak tape sk P= A, x 


表示 纯 基 体 梁 的 刚度 系数 。 计 算得 到 两 端 夹 紧 边 界 
条 件 下 , 当 GPL 的 重量 分 数 分 别 为 0.1% ALO. 3% 时 
梁 的 无 量 纲 届 曲 临界 载荷 的 数值 结果 ,并 将 该 临界 
载荷 与 文献 14] 中 使 用 微分 求 积 法 所 获得 的 相应 
结果 同时 列 入 表 1 中 。 由 表 1 可见, 用 本 人 研究 中 的 
解析 求解 方法 所 得 的 结果 与 文献 [ 14 ] 中 数值 方法 
所 求 的 结果 极其 接近 , 故 证 明了 本 研究 理论 推导 及 
解析 求解 过 程 的 正确 性 。 
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表 1 石墨 烯 增强 功能 梯度 梁 的 无 量 纲 屈曲 临界 载荷 比较 


Tab.1 Comparison of dimensionless critical buckling 


loads for the FG-GPLRC beams 


Zcp,/%  U-GPLRC  X-GPLRC 0-GPLRC 
本 研究 0.1 0.038 4 0.0410 0.033 2 
文献 [14 0.1 0.039 1 0.043 3 0.0349 
本 研究 0.3 0.0577 0.0703 0.044 1 
文献 [14 0.3 0.058 6 0.0709 0.045 8 
首先 画 出 了 石墨 烯 增强 功能 梯度 梁 的 前 四 阶 屈曲 
模 态 (图 2) ,由 图 2 可 见 该 复合 材料 梁 的 屈曲 构 型 与 均 


1.0 


匀 材 料 或 者 一 dd alaia 


图 2 石墨 烯 增强 功能 梯度 梁 的 前 4 阶 届 曲 模 态 
Fig.2 The first four order buckling modes 
of the FG-GPLRC beam 
其 次 分 析 FG-GPLRC 梁 的 临界 届 曲 载荷 及 其 影 
响 因 素 。 图 3 给 出 了 不 同 石墨 烯 重量 分 数 梁 的 临界 
屈曲 载荷 ,其 中 横 坐 标 和 纵 坐 标 分 别 表 示 GPL 的 宽 
FEE EAE KS Il AJE HH ER TT o 


400 


X-GPLRC 
350b Lon /Don=1, Uh=10 


Don /hen 
图 3 石 浴 燃 质 量 分 数 对 临界 届 曲 载荷 的 影响 


Fig.3 Effect of graphene mass fraction 


on the critical buckling loads 
由 图 3 可 见 当 GPL 在 基体 中 的 占 比 越 高 时 , 梁 
的 临界 屈曲 载荷 也 越 大 。 尤 其 当 石 墨 烯 片 的 宽厚 比 
较 大 时 , 增强 效果 更 加 明显 ,石墨 炳 片 宽厚 比 为 
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10° .含量 仅 为 0.3% 时 的 FG-GPLRC 梁 的 临界 届 曲 
载 向 是 纯 基 体 梁 的 4.2 倍 。 

为 了 更 进一步 考察 石墨 烯 片 的 几何 尺寸 对 梁 届 
曲 的 影响 ,图 4 画 出 了 梁 的 临界 届 曲 载荷 随 GPL 的 
几何 参数 变化 的 关系 曲线 ,其 中 图 4(a) 和 图 4 Cb) 
的 横 坐 标 分 别 表示 GPL 的 宽厚 比 和 长 宽 比 , 纵 坐 标 
为 梁 的 临界 届 曲 载荷 。 


400 


X-GPLRC 
Sam=0.3%, Up=10 


350 


-一 一 = 


Pl 
> 


250 
200 
= 
> 
150 
CD o 7 。 
< Dis. Mess: 
= (a) 宽厚 比 的 影响 
© 280 
© BonW03%, banhan =100 sn nnnnnnn= 
© 2% 
at i ee ee gaara] 
= ne A 
N 
E u 
a ---+X-GPLRC 
; = 一 = O-GPLRC 
>< 
200% 
二 o 5 T T a 
: = lonl bon 
= (b) 长 宽 比 的 影响 


Qaa 石墨 烯 片 的 几何 尺寸 对 临界 届 昌 载荷 的 影响 
Fig.4 The influence of the geometry size of 
CPL on the critical buckling loads 

在 石墨 烯 总 含量 和 石墨 烯 的 宽度 不 变 的 情况 
下 , 较 大 的 长 宽 比 意味 着 石墨 烯 片 具有 和 较 大 的 表面 
积 , 较 大 的 宽厚 比 意味 着 在 石墨 燃 中 所 含 的 单 层 石 
墨 烯 片 较 少 '% 。 因 此 ,从 图 4(a) 中 可 以 发 现 随 着 宽 
厚 比 的 增 大 ,临界 屈曲 载荷 也 随 之 增 大 ,说明 含有 较 
少 的 单 层 石墨 烯 片 时 增强 效果 更 加 明显 。 但 当 宽 厚 
比 大 于 10 时 这 种 变化 不 太 明 显 , 即 宽厚 比 大 于 10° 
时 ,再 无 法 通过 增 大 GPL 的 宽厚 比 来 有 效 提高 梁 的 
强度 。 从 图 4(b) 中 看 到 随 着 长 宽 比 的 增 大 ,临界 届 
曲 载荷 也 随 之 增 大 ,说 明 具 有 较 大 表面 积 的 石墨 烯 
片 对 于 功能 梯度 梁 的 强度 具有 明显 的 增强 效果 。 但 
当 lop/bem > 4 后 对 临界 屈曲 载荷 的 影响 变 得 不 太 


明显 , 即 lop / ben, > 4 时 ,也 无 法 通过 增 大 GPL 的 长 
宽 比 来 提高 增强 效果 。 由 图 4(b) 也 发 现 ,不 同 的 
GPL 分 布 形式 对 临界 屈曲 载荷 的 增强 效果 不 同 ,X- 
GPLRC 模式 能 更 有 效 地 提高 梁 的 临界 载 答 ,增强 效 
果 最 好 , 即 在 梁 的 上 下 表面 分 散 由 较 少 单 层 石 墨 烯 
片 组 成 的 石墨 烯 是 提高 粱 的 临界 载荷 的 最 有 效 
方法 。 

最 后 分 析 FG-GPLRC 梁 的 后 届 曲 响应 及 其 影响 
因素 。 图 5 和 图 6 分 别 绘 出 了 GPL 的 质量 分 数 和 分 
布 形式 不 同时 的 石墨 烯 增强 功能 梯度 梁 的 后 届 曲 路 
径 , 其 中 横 坐 标 和 纵 坐 标 分 别 表示 梁 的 轴 向 载荷 和 
后 届 曲 最 大 接 度 。 由 图 5 可 见 , 随 着 GPL 质量 分 数 
的 增 大 , 梁 的 后 届 曲 承载 能 力也 随 之 显著 增 大 。 由 
图 6 可 见 ,3 种 分 布 模式 下 的 FG-GPLRC 梁 都 具有 
稳定 的 后 届 曲 路 径 , 并 且 X-GPLRC 分 布 模式 下 梁 具 
有 更 高 的 后 届 曲 承载 能 力 , 其 次 是 U-GPLRC 分 布 ， 
最 后 是 O-GPLRC 分 布 。 也 就 是 说 X-GPLRC 分 布 模 
式 在 上 表面 和 下 表面 具有 更 高 的 GPL 含量 , 梁 的 弯 
曲 刚 度 最 高 ,所 以 梁 具有 更 强 的 承载 能 


4.0 


X-GPLRC, Wi=10 
lom/bom=1, Bon/han=1 000 


w(0.5) 
2 


1 
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图 5 A2 n BIE E OT E J a 8 1 


Fig.5 Effect of the graphene mass fraction 
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on the post-buckling paths 
25 
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K6 3 种 GPL 分 布 模式 梁 的 后 届 曲 路 径 
Fig.6 Post-buckling paths of the beams 
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with three GPL distribution modes 
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图 7 和 图 8 表示 了 GPL 的 几何 尺寸 对 后 屈曲 路 
径 的 影响 ,图 中 不 同 曲线 分 别 代 表 GPL 的 长 宽 比 和 
宽厚 比 不 同 梁 的 后 屈曲 路 径 。 由 图 可 见 , 随 着 GPL 


张 靖 华 ， 等 ; 功能 梯度 石墨 烯 增强 


nay iAd HE: 
afke ia t am VOE sos 


由 其 变化 曲线 可 知 不 同 GPL 4AE PR IN AY) Je JH 
曲 变形 有 一 定 的 差异 ,这 是 因为 不 同 分 布 模式 梁 的 
强度 不 同 , 相 比 其 他 分 布 模式 ,X-GPLRC 分 布 模式 


的 长 宽 比 和 宽厚 比 的 增 大 , 深 的 后 屈曲 最 大 挠 度 逐 
渐 减 小 ,具有 更 高 的 后 屈曲 承载 能 


4.0 


3.5 X-GPLRC 
Ban=0.3%, Bem /hon.=10 
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=== lon/bon=10 
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Fig.7 Effect of geometry size of graphene platelets 


on the post-buckling paths 


X-GPLRC 
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Fife) Effect of width-to-thickness ratio of graphene platelets 


naXiv 


on the post-buckling paths 
图 9 绘 出 了 不 同 分 布 形 式 的 石墨 烯 增强 功能 梯 
度 梁 的 后 届 曲 构 型 图 ,其 中 给 定 无 量 纲 载 寄 P=250。 


1.0 
— X-GPLRC 
a CPLRG e 
08 IEE O-GPLRC re s besi /ham=100 


图 9 FG-GPLRC 梁 的 后 届 曲 构 型 
Fig.9 Post-buckling configurations of the FG-GPLRC beams 
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本 研究 通过 精确 解析 求解 非 线 性 微分 方程 研究 
了 石墨 烯 增强 功能 梯度 梁 的 届 曲 及 后 届 曲 特性 , 获 
得 了 屈 曲 临界 载荷 `. 届 曲 模 态 及 后 屈曲 挠 度 的 解析 
表达 式 。 并 通过 解析 表达 式 的 数值 算 例 探讨 了 不 同 
的 GPL 分 布 模式 、 几 何 尺 寸 和 重量 分 数 等 参数 对 石 
墨 烯 增强 功能 梯度 梁 的 非 线 性 届 曲 及 后 届 曲 的 影 
响 , 获 得 如 下 主要 结论 。 

1) 加 入 极 少 量 的 GPL 可 以 显著 提高 功能 梯度 
梁 的 临界 屈曲 载荷 和 后 屈曲 承载 力 , 即 对 梁 的 强度 
和 刚度 具有 明显 的 增强 效应 。 

2) 当 GPL 具有 较 大 的 表面 积 和 较 少 的 单 层 石 
黑 烯 片 时 ,能 提高 对 梁 的 增强 效应 ,其 临界 载荷 增 
大 ,后 屈曲 挠 度 减 小 。 

3) 不 同 的 GPL 分 布 形 式 对 梁 的 增强 效果 不 同 ， 
当 石 黑 烯 以 X-GPLRC 分 布 时 , 即 梁 的 上 下 表面 具有 
较 多 的 GPL, 梁 具有 更 高 的 承载 能 力 以 及 临界 届 曲 
BAW o 
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